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THERMAL DISTORTIONS OF THIN-WALLED PARABOLOIDAL SHELLS* 


By Joseph E. Walz 
Langley Research Center 

SUMMARY 

Equations are presented to describe the linear behavior of a deep thin-walled shell 
of revolution subjected to an axisymmetric thermal loading. A solution is obtained for a 
paraboloid with a ring at the edge. Results are plotted showing the rotation of the tangent 
to the shell due to (1) constant temperature difference between shell and ring and (2) con- 
stant temperature gradient through the thickness of the shell. 

INTRODUCTION 

The radiation of the sun provides an attractive source of power for satellites and 
space vehicles. To transform the radiation of the sun into electrical power, a conversion 
system is required. In order to obtain usable amounts of electrical power, some con- 
version systems require concentration of the rays of the sun; thus, the concept of a solar 
concentrator has been developed. Several shapes of concentrators are conceivable, but 
most of the attention has been devoted to the paraboloidal shape because of its high con- 
centrating ability. In the utilization of a paraboloidal mirror as a solar concentrator, the 
question arises as to the effect that temperature gradients have on the shape of the shell 
forming this mirror and consequently on its concentrating ability. 

In order to obtain maximum efficiency, efforts will be made to train the axis of the 
paraboloidal mirror directly at the center of the sun; hence, axisymmetric temperature 
effects may be encountered. In preliminary studies (such as ref. 1) where the use of 
linear plate formulas have been employed, temperature gradients through the thickness 
of the paraboloid were shown to produce large overall distortions. Reference 2 develops 
a means of solving the general axisymmetric thermal problem for an isotropic shell of 
revolution which is based on deep shell theory; unfortunately, the numerical results in the 
illustrative example section appear to be in error. 

The information presented herein is based in part upon a thesis submitted in 
partial fulfillment of the requirements for the degree of Master of Science in Engineering 
Mechanics, Virginia Polytechnic Institute, Blacksburg, Virginia, June 1965. 



In the present paper equations governing the linear behavior of a general isotropic 
shell of revolution subjected to an axisymmetric temperature distribution are given. 

Deep shell theory, although different from that used in reference 2, is employed to obtain 
these equations. These equations are specialized for a paraboloid with a ring at the edge 
subjected to a constant overall temperature change and a constant temperature gradient 
through the thickness; the ring may have a uniform temperature change which is different 
from that of the shell. The ability of the ring to restrict the expansion and rotation of the 
edge of the shell can be varied by adjusting the stiffness parameters involved. Results 
are presented in equation form for stresses and displacements, and figures showing 
rotation of the tangent at any location are plotted for a wide range of parameters. An 
indication is given as to the effect that thermal distortions have on the concentrating 
ability of a paraboloid used as a solar concentrator. 

SYMBOLS 


U.S. 


The units used for the physical quantities defined in this paper are given both in 
Customary Units and in the International System of Units (SI). (See ref. 3.) 


Ar 

a 


b 

C 



D 

E 

E r 


cross-sectional area of ring 


radius from axis to centroid of ring (see fig. 3) 


cross-sectional radius of circular ring (see fig. 3) 


extensional stiffness of shell, Eh 


dimensionless ratio of extensional stiffness of ring to that of shell, 


E r A 


r^ r 


Ehr f 


concentration ratio herein used as ratio of projected area of concentrator to 
area of aperture of an absorber located in focal plane 


flexural stiffness of shell, 


Eh J 

12(1 - i/2) 


dimensionless ratio of flexural stiffness of ring to that of shell, — 

r 0 D 

Young’s modulus of shell 
Young’s modulus of ring 


2 



e h 


e v 



f 

H 

h 

Ir 



ijj 


Kl,K 2 

k 


l 

*a 

M|,M 0 

m 

N ?J N 0 


P 


radial eccentricity of centroid of ring from point of attachment of ring to 
shell (see fig. 3) 


axial eccentricity of centroid of ring from point of attachment of ring to 
shell (see fig. 3) 


elliptic integral of first kind, 

r o 

focal length of paraboloid, — 
radial stress resultant 



2 

o 


f 


d X 

- k^sin^y 


thickness of shell 


moment of inertia of cross section of ring about centroidal axis in plane 
of ring 


integers 

constants 

modulus of elliptic integral 

distance in focal plane that a reflected ray misses focal point 
radius of aperture of absorber 

bending moments per unit length of sections of shell (see fig. 1) 
applied uniform moment per unit length 

normal stress resultants in meridional and circumferential directions, 
respectively (see fig. 1) 

applied uniform force per unit length 


Ph 


radially distributed loading, positive in +r-direction 


p v axially distributed loading, positive in +z-direction 

Q transverse shear stress resultant 

r radial distance to point on meridian of shell 

r 0 radial distance to edge of shell 

T temperature change from stress-free datum 

AT = Tanner - T ou t er 

T^ner temperature change from datum of point on inner surface of shell 

Touter temperature change from datum of point on outer surface of shell 

T r unif orm temperature change of ring from datum 

m T ou t er + Tinner 

Ts = 2 

u radial displacement of shell 

u nondimensional radial displacement of shell 

u a total radial displacement of attachment point of ring 

u aci> u ac2 ra -dial displacements of point of attachment relative to centroid of ring 

u cl> u c2 radial displacements of centroid of ring 

V axial stress resultant 

w axial displacement 

w nondimensional axial displacement 
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z axial distance to point on meridian of shell 

z Q axial distance to edge of shell 

a Lame parameter defined by equation (13a); specialized for paraboloid by 

equation (38) 

Q! r coefficient of linear thermal expansion of ring material 

0 's coefficient of linear thermal expansion of shell material 

/3 rotation of tangent of point on middle surface of shell from undeformed 

position to deformed position (see fig. 1) 

y rim angle 

Kronecker delta (6jj = 0, if i ^ j; Sy = 1, if i = j) 

e^eg strain in meridional and circumferential directions, respectively, at any 

point in shell 

strain in meridional and circumferential directions, respectively, at points 
on middle surface of shell 

C coordinate along inward normal of shell (see fig. 1) 

t ? expression defined by equation (41); specialized for paraboloid by 

equation (42) 

9 circumferential coordinate 

curvature changes in meridional and circumferential directions, respectively 
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X 


parameter indicating amount of concavity of paraboloid, 



V 2 *. 


v Poisson's ratio of shell material 

| dimensionless coordinate specifying point on meridian, — 

r o 

p expression defined in equations (28); specialized for paraboloid by 

equation (39) 

normal stress in meridional and circumferential directions, respectively 

$ expression defined in equations (28); specialized for paraboloid by 

equation (40) 

<p angle that tangent to a point on undeformed middle surface of shell makes 

with plane perpendicular to axis of shell (see fig. 1) 

X amplitude of elliptic integral 

Subscript: 

o evaluated at £ = 1 

Primes denote derivatives with respect to £. 

LINEAR AXISYMMETRIC SHELL EQUATIONS 

The governing equations for a linear thermoelastic axisymmetric problem involving 
a shell of revolution are obtained by modifying the equations in reference 4 to include 
temperature effects. The positive directions of the stress resultants and couples, dis- 
placements, and rotation involved in the development of these equations are shown in 
figure 1. 


Stress-Strain Relations 

When transverse normal stress is neglected, the following stress strain relations 
hold as long as E, Q! s , and v are taken to be constant over the range of temperatures 
being considered: 
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( 1 ) 


(1 - )a ^ + ve^j - Eo! s T(1 + v) 

(1 - v2)o e = E^e e + ve - Ea s T(l + u) (2) 

The assumption is made that the temperature is linear through the thickness of the shell; 
thus 


r 

T = T s + — AT 
s h 

where 


(3) 


-K 


T inner + ^outer) 


( 4 ) 


AT - Tj nner - T ou ter (5) 

The following strain variations through the thickness are used (ref. 4) when attention is 
restricted to thin shells (that is, ratio of shell thickness to radius of curvature is small 
compared with one): 


e a- e z + 

e o = 6 e + % K eJ 


(6) 


In obtaining the relationships expressed in equations (6), the deformations due to trans- 
verse shear and transverse normal stress are neglected; thus, points on the normal to 
the undeformed surface lie on the normal to the deformed surface. The stress resultants 
and couples are defined as 


C h / 2 

N i'J-h/2 a « dC 


N, 


= f 


h/2 


On d? 


-h/2 
C h /2 
* ■ 3 _ h/2 

r h/2 

Mo = \ CToi: d 5 

6 -3 -h/2 9 




m 
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Substitutions from equations (1), (2), (3), and (6) into equations (7), after performing the 
required integrations, yield 


where 


xt C f— — \ CofeTs 
Nt = let + veA 

5 1 - v2\ % 6 ) l-i/ 


C /- - \ Ca s Ti 

N Q - f 6 a + ve fc) 

i - A e v i-v 


l-i/ 

Coj S Ts 
■ v 

+ vkq\ - D(1 + i/)-^-j- — 
M e = d(kq + VK^j - D(1 + v)— 




( 8 ) 


C = Eh 


D = 


Eh' ; 


( 9 ) 


12(1 - i/2) 




Equilibrium and Compatibility Equations 

Equilibrium of forces in the axial and radial directions is found in reference 4 
to be (H is the radial force resultant; V is the axial force resultant (see fig. 1)) ; 


(rV)' + ropy = 0 

(10) 

(rH)' - 3N q + r5p h = 0 

(ID 


and the moment equilibrium equation is 


^rM^' - oMg cos <p + ra(H sin cp - V cos (p) = 0 


where a prime denotes differentiation with respect to and 

1 1/2 


5= j(r’) 2 +(z’) 2 1 


tan cp - =— 

Y' 


( 12 ) 


(13a) 

(13b) 
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A compatibility relation is given which for small rotations j3 is 


(re^ - r T e ^ = r T /3 tan cp 

Strain Displacement Relations 
Strain displacement equations given in reference 4 are 

— u 


and for small rotations 


— w T + r ' B 
— z' 


( 14 ) 


(15) 


(16) 


where u and w are radial and axial displacements, respectively. Thus, substitution 
of strains from equations (8) into equations (15) and (16) yields 




+ r «s T s 




z'«s T s d £ 


The curvature changes are given as 





(17) 

(18) 


(19) 


and for small rotation 


B cos w 


( 20 ) 


Development of Governing Second-Order Differential Equations 

The complete system of field equations can be reduced to two ordinary second- 
order differential equations in two unknowns. To develop these equations, it is desirable, 
initially, to express all quantities in terms of H, V, /3, and the distributed and thermal 
loads. The normal stress resultant and the transverse shear stress resultant Q 

are related to the axial and radial stress resultants by the following equations: 


9 



N| = H cos cp + V sin cp 


( 21 ) 


Q = -H sin + V cos cp 

Once p v is established, the solution for V from equation (10) is 

V= -If r ^V d 4 


( 22 ) 


(23) 


Substitution of the moment expressions in equations (8) into the moment equilibrium 
equation (12) after first using the curvature changes given in equations (19) and (20) yields 
a second-order differential equation involving /3 and rH. The other differential equa- 
tion required is obtained by substitution of Tq and <T| from equations (8), with Ng 
and N| from equations (11) and (21), into the compatibility equation (14). By use of the 
following transformation, 



(24) 

(25) 


the two differential equations so derived become 


where 


X - 0X + 2 


% y * - * a * ifj /2 h^ 


Y" - ipY - 2f-^-jp2$ 2 X = G 2 - (r5C) 1/2 (a s T s y 


(26) 

(27) 


p 2 = 3(1 - v 2 ) 


z'a 


rh 


ref 


$ 


2 - 


z J a 


rh 


z 'a 

“rh 


ref 


j 


(28) 

Equation continued on next page 
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v J 

an d | | re j means the absolute value of a quantity evaluated at a suitably chosen refer- 
ence point on the shell. 


Approximate Solution to Governing Equations 

It is convenient to consider separately the homogeneous and particular solutions of 
the equations. In equations (26) and (27) the quantity p2 is, in general, of the order of 
the ratio of a radius of curvature to the thickness of the shell and is large compared with 
unity. The functions 0 and i// appearing in these equations are, in general, of the 
order of unity provided the shell is such that the thickness is slowly varying (it may not 
vary appreciably over a distance on the order of ^Rh, where h is the average thickness 


11 






and R is a representative radius of curvature) and provided the functions are not eval- 
uated in a shallow region of the shell. Under these circumstances, for an approximate 
solution to the homogeneous parts of equations (26) and (27), it would be permissible to 
omit the 0 and the ip terms and put the equations in the following form: 

(x H + 1Y h ) H - 2i JL. p 2*2^x H + iY H ) = 0 (29) 


where the subscript H refers to the homogeneous part. Reference 4 seeks to obtain 
an asymptotic solution to equation (29) by considering the following series: 


X H = e p ^(xi + X 2 p -1 + X 3 p- 2 + . . .) (30) 

Y H = e p/i( ^Y! + Y 2 p -1 + Y 3 p~2 + . . j (31) 


As in reference 4 the solution of equation (29), except for small terms of order 1/p, is 
obtained by retaining only the leading terms in the series in equations (30) and (31) to 
yield 


X H + iY H = * 


_ *- 1/2 


C^e 


(l+i)p^ $ d£ ^ -(1+i )p^ $ d£ 


(32) 


z ^ 

where Cj and C 2 are complex constants, - — — is taken to be 1, and 

I z ' I 

±(1 + i) V $ d£. The case where — - — = -1 requires changing i to -i 
J | z» | 

nents in equation (32). 


ju is 

in the expo- 


In the treatment of approximate particular solutions, it is to be observed that the 
right-hand sides of equations (26) and (27) consist of applied pressure loading terms and 
thermal loading terms. With the stipulation that the derivative of the pressure loading 
with respect to £ is not of order' of magnitude greater than the loading itself, conven- 
tional membrane theory has been shown to provide an adequate approximate particular 
solution consistent with retaining only the leading terms in the homogeneous solution 
(ref. 5). To obtain a particular solution (denoted by subscript P) for the thermal loading 
under a similar condition, the first two terms on the left-hand sides of equations (27) and 
(26) would be neglected, and, as a result, 


Xp = 


1 

2 p 2 # 2 




(33) 
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which correspond to those obtained in reference 6 for constant thickness. As far as 
stresses are concerned, the retention of these particular solutions appears to be incon- 
sistent with retaining only the leading terms in the homogeneous solution and stresses 
obtained in this manner are in question; however, meaningful values for the rotation of 
the tangent may be obtained away from the edge of the shell. 

Separation of equation (32) into real and imaginary parts and the addition of the 
particular solutions yields the complete solution for X and Y. Thus, the deformations 
may be evaluated once the geometry is specified. 

PARABOLOIDAL SHELL WITH EDGE RING 

The general equations developed in the preceding sections are now specialized to a 
paraboloid cf revolution of constant thickness with a ring at the edge. In particular, 
there is no loading other than a thermal loading due to temperature changes, T s of the 
shell and T r of the ring, and constant temperature gradient through the thickness of 
the shell AT. The presence of a ring is considered because a paraboloid used as a 
solar concentrator may have a support ring attached at its edge, but other boundary con- 
ditions such as a free or clamped edge may be realized by suitably adjusting two stiffness 
parameters of the ring. 

The dimensionless coordinate £ is used to specify a point on the meridian where 

£ = (35) 

r o 

and r is the radial distance from the axis to the particular point on the meridian and 
r Q is the radial distance to the edge. In this coordinate system a paraboloid is specified 
by 

z -^=! xl/2r °« 2 (36) 

where z is the axial location of a point, f is the focal length of the paraboloid, and 
. The rim angle y and X have the relationship expressed by the following 
equation, and a plot of this relationship is shown in figure 2: 




(37) 


tan y = 


2A 1 / 2 

1 - X 


The rim angle y is the angle formed between the axis of the paraboloid and a line 
passing through a point on the rim and the focus (fig. 2). 

The Lame parameter, given by equation (13a), for a paraboloid in this coordinate 
system becomes 


a = r 0 \[l + A£ 2 (38) 

and the first two quantities of equations (28) become 

p= 2) k ^ (39) 

$ = (1 + X| 2 ) 1//4 (40) 

when the apex is selected as the reference point. Another quantity r/ defined as 


V = P 



in this coordinate system, becomes 


V = P 






(41) 


(42) 


where F(k,y) 



is an elliptic integral of the first kind and 


X- 



(43) 


The particular solutions Xp and Yp, as given by equations (33) and (34) are both 
zero, since no applied loading is considered besides thermal loading, the shell is of con- 
stant thickness, and the temperature components are constant. Thus, the homogeneous 
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solutions of equations (26) and (27) are the complete solutions. For the case of a parab- 
oloid which is closed at the apex, C 2 is equal to zero (ref. 7) and the following solution 
is obtained: 


X = ^Eh 3 4? -1 / 2 e ^(Ki cos 77 - K 2 sin rij 

(44) 

Y = ^Eh 3 4> - * //2 e sin 77 + K 2 cos 77 j 

(45) 


where and K 2 are dimensionless constants obtained by nondimensionalizing the 

real and imaginary parts of Cj. Henceforth, as in equations (44) and (45), the sub- 
script o indicates a quantity evaluated at the point £ = 1. Equations (44) and (45), 
with the use of equations (24) and (25), permit the evaluation of stresses and deformations 
once the constants and K 2 have been evaluated from the boundary conditions. 

Deflection and Rotation of Edge of Shell 

Boundary conditions for a shell having a ring attached at its edge are such that 
deformation and force compatibility must be maintained at the attachment point of the ring 
and the shell. Thus the radial displacement u of the shell at the edge must be equal to 
the radial displacement of the ring at the point of attachment. Also, if the method of 
attachment permits transfer of moment, the edge of the shell and the ring must rotate the 
same amount about an axis perpendicular to the cross section of the ring. 

For the shell, from equations (24) and (44), the rotation at any point is 

(3 = 2^3(1 - 1 , 2 ) (l + x ^) 1 / 8 !- 1 / 2 ^ 770 cos 77 - K 2 sin 77 ) (46) 

and thus at the edge of the shell 

/3 0 = 2^3(l - i/2)(l + X) 1 / 8 ^ cos rj Q - K 2 sin (47) 

where the subscript o denotes evaluation at £ = 1. An approximate relation for the 
u-displacement, based on the consideration that the solution of the governing equations 
neglects terms of order 1/p compared with 1, is found from equations (17), (11), (25), 
and (45) to be 

u = hp (l + A£ 2 ) ^ £ 1//2 e ^|~Ki(sin 77 + cos 77 ) + K 2 (cos 77 - sin 77 )j + r 0 £agT s (48) 


15 


I 



A nondimen sionalized radial displacement u is obtained by dividing equation (48) by h; 
thus, at the edge of the shell 


u 0 = (1 + X)" 1/8 p K^sin r] 0 + cos r] Q "j + K 2 ^cos rj 0 - sin r) Q 




■)] + T 


“s T s 


(49) 


Deflection and Rotation of Ring 

For a ring subjected to a uniform radial loading per unit length P acting in the 
plane of the ring, there is no rotation; and the displacement u £l in the direction of 
loading is given by 


u 


cl 


Pa 2 

E r A r 


(50) 


(See ref. 8.) For a ring subjected to a uniform moment per circumferential length m 
tending to twist the ring inside out, there is no displacement but a rotation in the direc- 
tion of the moment such that 


/3r = 


ma 2 

E 


(51) 


(See ref. 8.) If the ring has attained a uniform temperature increase T r from initial 
equilibrium temperature, the centroid of the cross section of the ring has an additional 
displacement 

u c2 = (52) 

Since, in general, the ring is attached at some point other than its centroid, addi- 
tional effects must be considered. Because the radial force does not act necessarily 
through the centroid, the moment created by this eccentricity must be included in equa- 
tion (51). Rotation of the ring produces a relative displacement between the centroid and 
the point of attachment. For small positive angular rotation /3 r , the radial displacement 
of the point of attachment relative to the centroid is 

u acl = e v&r (53) 

where ey is positive when the point of attachment is above the centroid as shown in 
figure 3. The displacement of the point of attachment relative to the centroid due to T r 
is 
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u ac2 — e h a r T r 


(54) 


where e^ is positive when the centroid of the ring is at a distance greater than r Q 
from the axis of the paraboloid. Thus, the total radial displacement u a of the point of 
attachment is given by 


u a = u cl +u c2 +u acl +u ac2 < 55 > 

Evaluation of Constants 

In the process of evaluating the constants Ki and K 2 there are three independ- 
ent temperature terms to consider. It is possible, because of the linearity of the prob- 
lem, to consider each effect separately and then combine these effects to obtain the total 
effect. Thus, 


K 1 = K ll(“sT s ) + K 12 (a s AT) + K 13 (« r T r ) (56) 

K-2 = K2i^Q!sT s ) + K 2 2 ^Q!s AT 1 ) + K23^0!rT r ) (57) 

To insure compatibility of the shell-ring structure, the nondimensional radial dis- 
placement of the edge of the shell determined from equation (49) is set equal to the non- 
dimensional radial displacement of the point of attachment of the ring determined from 
the proper nondimensionalization of equation (55). Similarly, the rotation at the edge of 
the shell determined from equation (47) is set equal to the rotation of the ring determined 
from equation (51). From force compatibility considerations, the quantities P and m 
appearing in equations (50) and (51) are 

Pa = -H 0 r 0 (58) 

ma — — ^4^ + HqCy^ r^ (59) 

(See fig. 3.) After the appropriate quantities from equations (8), (19), (25), (45), (46), 

(50), (52), (53), (54), (58), and (59) have been substituted into equations (55) and (51), the 
two aforementioned equalities result in two simultaneous equations which have the form 

gll K lj + g 12 K 2j = Ilj (60) 
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g 21 K lj + g 22 K 2j = g 2j 


(61) 


where the subscript j takes on the value 1, 2, or 3, depending on the temperature com- 
ponent of interest and 

~\ 

g ll = ~( r°r o eh ) sin + Y 2 ! 3 ! 1 ’ u2 ) cos Vo - PCI 1 + A) -1 / 4 (sin + cos 


g 12 = ~( ~Y~ ) cos V 0 “ c ^ Y 2 ^ 3 ( 1 " ^ 2 )sin r, Q - pC(l + X) _1 / 4 (cos r] 0 - sin 


Vo 


g 2i = p(sin t) q - cos 77 0 )(1 + ^)' 1/4 " ^ Y 2 \ 3 ( X “ ^ Sin V ° ' D (? ~ Te h ) cos v o 
g 22 = p(sin Vo + cos 77 0 )(1 + ^)" 1/4 - ^ Y 2 \( 3 ( 1 " " 2 ) cos V o + 

Sjj = C(1 + X) ~ QJgTg 


sin r) 0 


g 12 


= 0 


g 13 = -C(l +X) _1/8 ^a r T r 


g 21 


= 0 


«22 - ~ T* 1 + -) (1 . + X) ~ 1/8 «s ^ 


f(i - * 2 ) 


S23 


= 0 


(62) 


where 


— E r A~ 
C = -£-J£- 
Ehr n 


- E_I 


D = 


r x r 

r 0 D 


(63) 
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Within the limits of the accuracy of the development and solution of equations (26) 
and (27), the deflections, stress resultants, and couples become approximately 


Eh - " (t) + X ^ 2 ) / r 3/2 e (r? ° ^(Kij sin 77 + K 2 j cos 77) 

= (l + X£ 2 ) ^ X 1 / 2 4 " 1//2 e ^ n ° ^(Kij sin 77 + K 2 j cos 77) 

|Kij(sin 77 


Q. 

Eh 


4 -l/2 e -(7 0 -^) 


+ cos 77 ) 


+ K 2 j(cos 77 - sin 77)J 

= |(t) K 1 " ^ 2 )] a - 1/4 ( 1 + M 2 ) ^ r 1/2 e (??0 ^[Kij(cos 77 - sin 77) 


- K 2 j(sin 77 + cos 77 )J - ^ 


AT 


6 2j 


12(1 - v) 

rV2 - --l/4 1/4/ . « \“l/8 


w 


(64) 

(65) 


( 66 ) 


(67) 


( 68 ) 


[ 3 ^ " ^ 1 ^ 4 (l + *4 2 ) r 1/2 e ^jKij(cos 77 - sin 77) 

, 1 a s AT 

- K 2 j(sin 7, + cos 77 )J - — _ — 6 2j 

£ = u, = i^j /2 [ 3 (l - , 2 )J 1/4 x 1 / 4 (! + ^ 2 )- 1 / 8 4 l/ 2 e -(^o^) [ Klj(sin „ + cos v) 

+ K 2 j(cos 77 - sin 77)j + ^ |a s T s 6 ij 
£ = w, = -(^) 1/2 [ 3 (l - ^)] 1/ V 1 / 4 ( 1 + x^2) 1/8 r l/2 e -(V^)[K lj( cos 77 + sin 77) 

+ K 2 j (cos 77 - sin 77)] + | ^ A 1//2 | 2 o!sT s 6ij + ( 70 ) 

/ 3 j = 2 ^ 3^1 - ^ 2 )(l + ^ Olo ^ cos 77 - K 2 j sin 77^ ( 71 ) 


(69) 
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where 6ij, called Kronecker delta, is defined as 

6 ij = 1 (i = j)| 

> (72) 

6 ij = 0 (i / j)J 

and constants Wj representing nondimensional rigid body translations are to be deter- 
mined from the boundary condition on w. The total value of a quantity is therefore the 
sum of the components (for example, /3 = + 02 + £ 3 ). 

RESULTS AND DISCUSSION 

Stresses and deformations throughout a paraboloid shell with a ring attached at the 
edge can be calculated from equations (64) to (71). Because the rotation of the tangent /3 
has the most significant effect on solar concentrator performance, results are presented 
for the variation of /3 with radial distance for two types of temperature distributions 
and various values for the structural and geometric parameters. The temperature dis- 
tributions represented are ( 1 ) uniform temperature change of the shell and a possibly 
different temperature change of the ring from the datum temperature (consideration of 
a r T r - a s T s ), and (2) uniform temperature gradient through the thickness of the shell 
with its middle surface and the ring at the datum temperature (consideration of Q! s AT). 
Finally, results are presented of simplified efficiency calculations described in the 
appendix to give some idea of the influence of these thermal loadings on the performance 
of a solar concentrator. 

Effect of Temperature Difference Between Shell and Ring 

The rotation of the tangent to the shell /3 as a function of radial distance r/r 0 is 
shown in figures 4 to 9 for temperature distribution o! r Tr " ^sTs- All results are for 
X = 1/3 (except fig. 7), r 0 /h = 5000 (except fig. 8 ), and no ring eccentricities (except 
fig. 9). In figures 4 to 9 for the case of the same material in the ring and the shell and 
the same uniform temperature change in the ring and the shell, no rotations occur any- 
where. It must be borne in mind that the slope changes presented in this set of figures 
are those of a material point of the shell and are not necessarily those of a spatial coor- 
dinate. In temperature distribution a! r T r - a s T s the T s component causes the shell 
to experience an expansion, and at a spatial coordinate there is a slope change because 
one material point is replaced by another material point which has a different slope than 
the original point. In this instance most of the paraboloid has expanded to the shape of a 
new paraboloid which has its focal length 1 + OgTg times as large as the original focal 
length. 
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Results for various values of the relative extensional and flexural stiffness param- 
eters C and D for a constant ratio of D to C (d/C = 10 000) are shown in figure 4. 
For small values of C and D the ring stiffnesses are small relative to those of the 
shell, and the shell is relatively unaffected by the presence of the ring. As the values of 
C and D are increased, the region near the edge of the shell is beginning to be affected 
by the presence of the ring. When the values of C and D approach infinite values, 
the edge of the shell becomes "clamped." In this condition rotation of the tangent to the 
shell at the edge is not allowed, and the radial displacement is constrained to be the same 
as that of the ring caused by the temperature change of the ring. 

Results for a large value of the relative flexural stiffness parameter, D = 10 000, 
coupled with various values of C, and results for a moderately large value of the rela- 
tive extensional parameter, C = 0.1, coupled with various values of D are shown in 
figures 5 and 6, respectively. Effects of variations in the shape parameter X and var- 
iations in the radius-thickness ratio r 0 /h are shown in figures 7 and 8, respectively, 
for infinitely large values of both C and D. As the shell becomes flatter ( X 
decreasing or rim angle decreasing), the rotations become larger in magnitude every- 
where in the edge region and the edge effects are propagated farther into the shell 
(fig. 7). Thus from the point of view of keeping shell rotations small, larger rim angles 
are desirable. Finally, the influence of small axial eccentricities of the ring are illus- 
trated in figure 9. Effects of small radial eccentricities of the ring are not shown 
because they are negligible. 

Effect of Temperature Gradient Through Thickness of Shell 

The rotation of the tangent to the shell j3 as a function of radial distance r/r Q is 
shown in figures 10 to 12 for temperature distribution a s AT. The effects of various 
values of the relative stiffness parameters C and D for a constant ratio of 
d/c = 10 000 are shown in figure 10 for a shell for which X = 1/3 and r 0 /h = 5000. 

As the values of both C and D approach infinitely large values, (3 approaches zero 
throughout the shell. The free edge case where C = D = 0, on the other hand, leads to 
the largest rotations. Effects of variation in the shape parameter X are shown in fig- 
ure 11 for r Q /h = 5000 for the free edge case. As in the case of temperature distribu- 
tion Oj.Tr - o s AT, larger values of X lead to reduced magnitudes of rotation. The 
results of linear flat-plate theory used for preliminary studies such as reference 1 for 
which 


r o r 
Og AT h r Q 


(73) 
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are not given in this figure since values obtained from this formula are found to be much 
larger and affect the entire shell. Finally, effects of variations in r Q /h are shown in 
figure 12 for A = 1/3 for the free edge case. 


Influence of Thermal Effects on Concentrator Performance 


The loss in efficiency due to thermal distortions based on the assumptions and 
analysis contained in the appendix is shown in figures 13 and 14. It is emphasized that 
these results are based on the assumption that solar energy is coming from a point 
source. Furthermore, interaction of thermal distortions with other features which 
reduce efficiency has not been incorporated. The effect on efficiency of temperature 
distribution of r T r - ci! s T s is shown in figure 13 when AT is zero, and the effect of 
temperature distribution Og AT when Q! r T r - a; s T s is zero is shown in figure 14. 


The effect of temperature distribution o r Tr - Q! S Ts on efficiency is shown in fig- 
ure 13 for a shell having values of A = 1/3, r Q /h = 5000, and v = 0.3. The shell has an 
attached ring with properties such that C = 1 and D = 106 which simulates one of the 
most detrimental conditions. The efficiency is shown as a function of the area concen- 
tration ratio, herein used as the ratio of the projected area of the paraboloid to the area 
of the aperture of an absorber in the focal plane, for values of |a! r T r - a s T s | equal to 


100, 50, 20, and 10 /iin./in. (jum/m). 


When a r ~ a s = 10 juin./in-°F 
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j um 
m-°K 


which 


is roughly the value of this coefficient for metals, these curves correspond to 10° F 
(5.6° K), 5° F (2.8° K), 2° F (1.1° K), and 1° F (0.56° K). 


The effect of temperature distribution a s AT on efficiency is shown in figure 14 
for a shell having values of r 0 /h = 5000, v = 0.3, and A = 1/3. The shell does not have 
an attached ring and this case leads to the greatest losses in efficiency due to a tempera- 
ture gradient through the thickness of the shell. Curves are presented for four values of 
o!g AT. Since the area in which thermal distortions occur is small, it is not surprising 
that the loss in efficiency is not great. 


CONCLUDING REMARKS 


Linear equations are presented to describe the action of a paraboloid under axisym- 
metric temperature loading consisting of a constant temperature change of the shell, con- 
stant temperature gradient through the thickness of the shell, and a constant temperature 
change of an attached ring. Slope distortions are found to be confined to a region near the 
edge of the shell and it is found that larger rim angles are desirable since they restrict 
the distortions to smaller magnitudes and to a smaller region. 
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For the case of a paraboloidal shell with a free edge subjected to a constant tem- 
perature gradient through its thickness, the results for preliminary purposes based on 
flat-plate formulas indicate the entire surface is affected and the rotations involved are 
much greater than those determined by the present analysis. It is found that an attached 
ring relieves the distortions due to a temperature gradient through the thickness of the 
shell. The attached ring however becomes the source of loading on the shell when the 
shell is at one temperature and the ring at another temperature. The larger the ratio 
of the extensional stiffness of the ring compared with that of the shell, the greater the 
distortions. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 19, 1966. 
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APPENDIX 


EFFICIENCY CONSIDERATIONS 

The objective of this appendix is to estimate the loss in concentrating efficiency due 
to thermal distortions. The concentrating efficiency used herein is defined to be the 
ratio of the specularly reflected energy directed within a prescribed area in the focal 
plane divided by the total specularly reflected energy (shadowing effects of the absorber 
being neglected). A point source of energy is treated. From this simplified analysis an 
idea might be gained as to the influence of temperature differences on the performance of 
a solar concentrator. The interactions of thermal distortions with other features which 
also reduce the efficiency of the collecting system such as, the sun’s finite image, 
irregularities in the concentrator surface, and so forth, have not been included. 

Consider the reflecting surface being distorted at a point such that a ray which 
originally was directed to the focus is now displaced a distance l in the focal plane. 

The point under consideration is shown in figure 15 after the reflecting surface has 
undergone a positive rotation. Attention here is restricted to a point within the region 
near the edge of the shell for which consideration of the effect of the rotation alone is 
sufficient since it can be shown that the effects of vertical and radial displacements are 
of higher order. 

When transverse shearing deformations are neglected as was done in the shell 
development, the rotation of the tangent to the reflecting surface is the same as the rota- 
tion of the tangent to the middle surface /3 which can be determined from equation (71). 
From the geometry of figure 15, 


l = r - (f - z)tan(2<p - 2/3) (Al) 

Nondimensionalization of equation (Al) by dividing by r Q yields 


— = /3 

r 0 sin 4 cp 


(A2) 


where geometrical considerations have been employed and in the expansion of 

tan(2 cp - 2/3) in powers of /3, /3^ and higher terms have been neglected. Equation (A2) 

may be rewritten as 


±= x ' 1/2 +X ^ 2 )' 

r °~ 1-X« a 


(A3) 
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where sin 4 cp has been expressed as a function of £, after cp has been determined 
through the use of equations (35), (36), and equation (13b) to be 


v = tan 


(A4) 


If an absorber has a flat circular aperture of radius Z a in the focal plane. 


^a _ _1 


(A 5) 


where is the area concentration ratio herein defined as the ratio of the projected 

area of the paraboloid to the area of the aperture of the absorber. At a particular point, 
r/ro ? if the absolute value of l / r G determined from equation (A3) is greater than 
l//C A , the reflected energy from that point is considered lost and if the absolute value 
of l /r Q is less than 1 /|^Ca, the reflected energy from that point is considered to be 
collected. For a particular value of the concentration ratio, the projected area of the 
concentrator from which energy is collected can be found graphically and the ratio of 
this area to the total concentrator projected area is the efficiency. The results shown 
in figures 13 and 14 are obtained by selecting several values of the concentration ratio 
and calculating the corresponding efficiency. 
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(a) Geometry of shell-ring structure. (b) Internal reactions at attachment point of ring to shell. 

Figure 3.- Geometry and internal reactions on an axial section of shell-ring structure. 
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Figure 8.- Effect of r 0 /h on rotation of tangent to shell near edge due to a temperature difference between ring and shell. AT = 0; constant T r ; constant T S ; v = 0.3; 

x = 1/3; C = D = « 
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Figure 13.- Effect of thermal distortions near edge of a paraboloid due to a temperature difference between ring and shell on concentrating efficiency as a function of 
concentration ratio. Source of energy is treated as a point source. AT = 0; constant T r ; constant Ts; v = 0.3; A = 1/3; r 0 /h = 5000; 0 = 1; D = 10 6 ; 

ey = 0 ; eh = 0. 
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Figure 14.- Effect of thermal distortions near edge of a paraboloid due to a temperature gradient through thickness of shell on the concentrating efficiency_as a function 
of concentration ratio. Source of energy is treated as a point source. Constant AT; a r T r - a s T s = 0; v = 0.3; A = 1/3; r 0 /h = 5000; C = 0; D = 0. 




Figure 15.- Effect of rotation of tangent to shell on path of a reflected ray. 
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